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(g) Dynamic bandwidth estimation and adaptation for packet communications networks. 



(57) Access control for a packet communications 
network includes a dynamic bandwidth updat- 
ing mechanism which continuously monitors 
the mean bit rate of the signal source and the 
loss probability of the connection. These values 
are filtered to remove noise and then used to 
test whether the values fall within a pre-defined 
acceptable adaptation region in the mean bit 
rate, loss probability plane. Values falling out- 
side of this region trigger bandwidth updating 
procedures which, in turn, result in acquiring a 
new connection bandwidth, and determining 
new fitter parameters and new parameters for a 
leaky bucket access mechanism. 
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Technical Field 

This invention relates to traffic management in packet communications networks and, more particularly, 
to traffic monitoring, traffic measurement filtering, and adaptive bandwidth adjustment for such networks. 

5 

Background of the Invention 

In order to avoid congestion and insure adequate traffic flow in packet communication networks, it is com- 
mon to control the access of packet sources to the network on an ongoing baste. I n order to successfully control 

10 traffic access, it is necessary, first, to accurately charac^rize the trafTic so as to provide appropriate bandwidth 
for carrying that traffic. Simple measurements which provide accurate estimates of the bandwidth require- 
ments of a source are taught in the copending appl ication Serial Number 07/942,873, filed September 10,1 992, 
and assigned to applicants' assignee. In this application, the parameters used to characterize traffic include 
R, the peak bit rate of the incoming traffic in bits per second, m, the mean bit rate of the incoming traffic in 

is bits per second, and b, the mean burst length of the traffic in bits. Rather than using the actual burst length, 
however, a so-called "exponential substitution" technique is used to calculate and equivalent burst length which 
would produce the same packet loss probability if the traffic were a well behaved exponentially distributed 
on/off process. For traffic widely differing from such an exponential process, this equivalent burst length pro- 
duces a much more accurate characterization of the actual traffic and therefore permits a higher density of 

20 traffic on the same transmission facilities. 

The measured parameters are used to control the access of signal sources to the network when the actual 
traffic behavior departs significantly from the initial assumptions. A leaky bucket mechanism is one technique 
for controlling access to the network when the traffic exceeds the initial assumptions, but yet permits trans- 
parent access to the network when the traffic remains within these initial assumptions. One such leaky bucket 

25 mechanism is shown in the copending application Serial Number 07/943,097, filed September 10, 1992, and 
assigned to applicant's assignee. More particularly, the leaky bucket mechanism of this application prevents 
saturation of the network by low priority packets by limiting the number of low priority packets which can be 
transmitted in a fixed period of time while imposing a minimum on the number of red packets transmitted at 
a given time. Such leaky bucket control mechanisms optimize the low priority throughput of the packet network. 

30 High priority-traffic, of course, is transmitted with little or no delay in the leaky bucket mechanism. 

The above-described mechanisms are suitable for controlling traffic only if the traffic is reasonably well- 
behaved and remains within the general vicinity of the initially assumed traffic parameters. The traffic man- 
agement system, however, must be structured to deal with traffic which is not well-behaved and which departs 
substantially from the initially assumed traffic parameters. If such a departure persists for any significant 

35 length of time, a new connection bandwidth must be assigned to the connection to accommodate the new traf- 
fic parameters. Such adaptation of the control system to radical changes in traffic behavior presents the prob- 
lems of filtering the traffic measurements to separate transient changes of traffic behavior from longer term 
changes, and determining reasonable ranges within which the initially assumed traffic parameters can be 
maintained and outside of which new connection bandwidths must be requested. A bandwidth too large for the 

40 actual traffic is wasteful of connection resources while a bandwidth too small results in excessive packet loss. 
Ancillary problems Include reasonable ease in implementation of the adaptation process and reasonable com- 
putational requirements in realizing the implementation. 

Summary of the Invention 

45 

In accordance with the illustrative embodiment of the present invention, dynamic adaptation of a traffic 
control system to changes in the traffic parameters is provided by defining a region within which adaptation 
is not required and outside of which a new bandwidth allocation must be requested. In particular, the bandwidth 
requirement is adjusted upward if the measurements indicate that either a desired maximum packet loss prob- 

so ability will be exceeded or if the traffic on that connection will start to unfairly interfere with other connections 
sharing the transmission facilities. The bandwidth requirement can be adjusted downward, on the other hand, 
if significant bandwidth savings can be realized for both the user of the connection and for the balance of the 
network, without violating any quality of service guarantees for all of the connections. In further accord with 
the present invention, these limits on the adaptation region are converted to values of effective mean burst 

55 length and mean bit rates. The measured effective mean burst length and mean bit rates are then filtered to 
insure that the filtered values are statistically reliable, i.e., that a sufficient number of raw measurements are 
involved to insure a preselected confidence level in the results. This minimum number of raw measurements, 
in turn, determines the amount of time required to collect the raw measurements, given the mean bit rate of 
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the traffic. This measurement time can be used to measure not only the statistics of the Incoming data stream 
to the leaky bucket, but also the effect of the leaky bucket on the incoming traffic. This latter measurement 
allows a measure of how well the leaky bucket is dealing with variances in the offered traffic and hence the 
packet loss probability. 

5 When the traffic parameters fall outside of the desired adaptation region, a new connection with a different 
bandwidth is requested in order to accommodate the changes in the traffic parameters. Compulations to de- 
termine the need for adaptation can be minimized by using extrapolation techniques between computed values 
marking the bounds of the adaptation region, and confining consideration to only the upper right quadrant of 
the adaptation region. 

10 The adaptation mechanism of the present Invention has the distinct advantage of insuring a continuously 
reasonable traffic management strategy by insuring dynamic adaptation when the needs of the connection or 
the needs of the network call for such adaptation. Moreover, unnecessary adaptation is avoided, reducing the 
overhead involved in such adaptation. 

is Brief Description of the Drawings 

A complete understanding of the present invention may be gained by considering the following detailed 
description in conjunction with the accompanying drawings, in which: 

FIG. 1 shows a general block diagram of a packet communications network in which the dynamic traffic 
20 management mechanism of the present invention might find use; 

FIG. 2 shows a graphical representation of a connection request message which might be used to set up 
initial connections and dynamically altered connections in the packet communications network erf FIG. 1, 
using the dynamic traffic management mechanism of the present invention; 

FIG. 3 shows a general block diagram of the network access subsystem for the network of FIG. 1 in which 
25 the dynamic traffic management mechanism of the present invention is implemented; 

FIG. 4 shows a block diagram of the traffic characteristics estimation and adaptation module forming a 
part of the network access subsystem of FIG. 3; 

FIG. 5 shows a graphical representation, in the mean bit rate-effective burst length plane, of the adaptation 
region outside of which new connection parameters are requested foran existing connection in accordance 

30 with the present invention ; 

FIG. 6 shows the adaptation region of FIG. 5 replotted on the mean bit rate, red marking probability plane 
to reduce the computational requirement of determining operation outside of this adaptation region; 
FIG. 7 shows a flow chart of the process for controlling the access to a packet communications network 
such as that shown in FIG. 1, using the adaptation region illustrated graphically in FIG. 6; and 

35 FIG. 8 shows a flow chart of the process for updating a connection when required by the procedure of the 
flow chart of FIG. 7. 

To facilitate reader understanding, Identical reference numerals are used to designate elements common 
to the figures. 

40 Detailed Description 

Referring more particularly to FIG. 1, there Is shown a general block diagram of a packet transmission sys- 
tem 1 0 comprising eight network nodes 1 1 numbered 1 through 8. Each of network nodes 11 is linked to others 
of the network nodes 11 by one or more communication links A through L Each such communication link may 

46 be either a permanent connection or a selectively enabled (dial-up) connection. Any or all of network nodes 
1 1 may be attached to end nodes, network node 2 being shown as attached to end nodes 1 , 2 and 3, network 
node 7 being shown as attached to end nodes 4, 5 and 6, and network node 8 being shown as attached to 
end nodes 7, 8 and 9. Network nodes 11 each comprise a data processing system which provides data com- 
munications services to all connected nodes, network nodes and end nodes, as well as providing decision 

so points within the node. The network nodes 11 each comprise one or more decision points within the node, at 
which point incoming data packets are selectively routed on one or more of the outgoing communication links 
terminated within that node or at another node. Such routing decisions are made in response to information 
in the header of the data packet The network node also provides ancillary services such as the calculation 
of new routes or paths between terminal nodes, the provision of access control to packets entering the network 

ss at that node, and the provision of directory services and topology database maintenance at that node. 

Each of end nodes 1 2 comprises either a source of digital data to be transmitted to another end node, a 
utilization device for consuming digital data received from another end node, or both. Users of the packet com- 
munications network 10 of FIG. 1 utilize an end node device 12 connected to the local network node 11 for 
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access to the packet network 10. The local network node 11 translates the user's data into packets formatted 
appropriately for transmission on the packet network of FIG. 1 and generates the header which Is used to route 
the packets through the network 10. 

In order to transmit packets on the network of FIG. 1, It is necessary to calculate a feasible path or route 

5 through the network from the source node to the destination node for the transmission of such packets. To 
avoid overload on any of the links on this route, the route is calculated in accordance with an algorithm that 
insures that adequate bandwidth is available for the new connection. One such algorithm is disclosed in the 
copending application, Serial Number 07/674,917, filed April 28, 1992, and assigned to applicant's assignee. 
Once such a route is calculated, a connection request message is launched on the network, following the com- 

10 puted route and updating the bandwidth occupancy of each link along the route to reflect the new connection. 
One such connection request message is shown in FIG. 2. 

In FIG. 2 there is shown a graphical representation of a connection request message to be launched from 
a source node in the network of FIG. 1 to a destination node in the network along a precalculated route. The 
connection message of FIG. 2 comprises a routing field 20 which includes the information necessary to trans- 
is mit the connection message along the precalculated route. Also included in the connection request message 
of FIG. 2 is a connection request vector 22 which characterizes the important statistical characteristics of the 
new packet source and which allows this new source to be statistically multiplexed with the previously existing 
signals on each link of the route. As wiO be discussed in detail hereinafter, the connection request vector in- 
cludes a relatively few parameters necessary to adequately characterize the packet source. As described in 

20 the afore-mentioned copending application, Serial Number 943,097 and assigned to applicant's assignee, 
these parameters might include the maximum bit rate for the source, the mean of that bit rate, and the equivalent 
burst length of packets from that source. 

The values in the connection request vector are used to test each link of the route to determine if the new 
connection can actually be supported by that link, and to update, separately for each link, the link occupancy 

25 metric for that link to reflect the addition of the new connection. If the link occupancy has changed since the 
route was calculated, the connection may be rejected at any node along the route, and the source node notified 
of the rejection. Finally, in FIG. 2, the control fields 23 include additional information used in establishing the 
connection, but which is not pertinent to the present invention and will not be further discussed here. Note that, 
when a connection is to be taken down, a connection removal message having the same format as FIG. 2 is 

30 transmitted along the route of the connection to be removed. The link occupancy of each link is then updated 
to reflect the removal of this connection by subtracting the metrics for the removed connection. 

In FIG. 3 there is shown a general block diagram of the source bandwidth management subsystem in ac- 
cordance with the present invention comprising a leaky bucket module 34 to which the user traffic on input 
line 39 is applied. The output of leaky bucket module 34 is applied to the network 10 of FIG. 1. A subsystem 

35 similar to that shown in FIG. 3 is provided for each source of user traffic to be applied to the network 1 0. These 
bandwidth management subsystems are located in the endnodes 12 of FIG. 1 and one such a bandwidth man- 
agement subsystem is provided for each direction of transmission between two communicating users. Al- 
though such bandwidth management subsystems can be realized with hard-wired circuit components, the pre- 
ferred embodiment utilizes a programmed computer since such an implementation is more readily modified 

40 to accommodate improvements and to reflect changes in traffic patterns. 

In leaky bucket access control module 34, packets are launched into the network on line 40 with one of 
at least two different priority classes, conventionally called "red" and "green," where green is the higher priority. 
Green packets are guaranteed a pre-specified grade of service based on an acceptable level of delay and loss 
probability within the network. The red packets do not have the same guarantees and are discarded before 

46 the green packets when congestion occurs. Bandwidth-conserving strategies for optimally marking packets 
in a leaky bucket mechanism is disclosed In the copending application Serial Number 943,097, filed September 
10, 1992, and assigned to applicants' assignee. In general, the function of the leaky bucket module 34 is to 
"shape" the traffic before it enters network 1 0 (FIG. 1 ), especially for user packet s not conforming to the initially 
provided statistical description, by marking such packets red. If the traffic characteristics stay within the initially 

so negotiated values, however, the red marking mechanism is sufficiently limited to Insure the promised loss prob- 
ability. If the incoming traffic characteristics depart substantially from the negotiated values, estimation and 
adaptation module 33 is invoked to take corrective action since the leaky bucket module 34 can no longer cope 
with the new traffic. 

As described in connection with FIG. 1, when a new connection is to be set up through network 10, an 
55 initial estimate of the traffic characteristics is made by the packet source. This estimate arrives at the band- 
width management system of FIG. 2 on line 36 together with the quality of service requirements on line 36. 
Such quality of service (QOS) requirements include acceptable loss probabilities, acceptable delays, real-time 
delivery requirements, and so forth. Connection agent 32 passes these connection requirements on to path 
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selection controller 30 which uses these requirements, together with the up-to-date network description in to- 
pology database 31 to calculate a connection path through network 10 (FIG. 1) which satisfies aD of these re- 
quirements. One optimum connection path selection controller is described in the copending application Serial 
Number 07/874,91 7, filed April 28, 1992 and assigned to applicants' assignee. Once calculated, the proposed 

5 connection path Is encoded In a connection request message such as that shown in FIG. 2 and launched as 
a bandwidth request onto the network 10 over line 37 of FIG. 3. 

The bandwidth request message of FIG. 2 traverses the calculated connection path and, at each node 
along the route, is used to reserve, in the next leg of the connection, the bandwidth required to satisfy the re- 
quirements of the connection request If sufficient bandwidth is available In each link of the connection along 

10 the computed path, the destination endnode 12 (FIG. 1) receives the request and transmits back an acceptance 
of the new connection. If. at any link along the route, insufficient bandwidth is available due to changes in the 
traffic patterns, a denial of the connection request is transmitted back to the source endnode. These bandwidth 
replies, whether negative or positive, are delivered back to connection agent 32 on line 38. If the connection 
is denied, the user source Is notified and another attempt at the connection can be made later. If the connection 

is is accepted, leaky bucket module 34 is activated and supplied with the appropriate parameters to control the 
access of the user traffic. The user then begins introducing traffic on line 39. At the same time, estimation 
and adaptation module 33 begins monitoring this incoming traffic to determine if any significant changes in 
the incoming traffic characteristics have occurred during the life of the connection. If so, module 33 notifies 
connection agent 32 to request a new bandwidth allocation, supplying connection agent 32 with the new traffic 

20 parameters required for the connection. As before, connection agent 32 launches a new bandwidth request 
on line 37 requesting the adjustment of the bandwidth of the connection. If the adjustment is accepted, the 
leaky bucket parameters are updated with the new traffic characteristics and estimation and adaptation module 
33 continues to monitor the incoming traffic, but with the new characteristics. Note that only a new bandwidth 
aDocation is requested, rather than a new connection. This saves the overhead involved in taking down the 

25 old connection and setting up a new connection. If the requested additional bandwidth is not available, the 
connection can be either taken down or given a lower priority, depending on the original negotiations with the 
sending party at the source node. The present invention is an improvement in the estimation and adaptation 
module 33 of FIG. 3. 

In FIG. 4 there is shown a block diagram of an estimation and adaptation module in accordance with the 
30 present invention comprising a box 41 which uses the initial traffic parameters, supplied on line 36 of FIG. 3, 
to calculate an adaptation region. For convenience, this adaptation region is measured in units corresponding 
to the initially assumed mean bit rate of the incoming user traffic and the initially calculated red marking prob- 
ability of the leaky bucket module 34 of FIG. 3. The actual input traffic is applied to mean bit rate filter 45 to 
determine the actual current mean bit rate of the incoming traffic. At the same time, the actual red marking 
35 probability of the leaky bucket module 34 (responsive to the actual traffic served by the leaky bucket module 
34) is applied to red probability filter 47. The function of filters 45 and 47 is to filter out the transient changes 
in mean bit rate and red marking probability. The filtered mean bit rate and the filtered red marking probability 
are compared, in unit 42, to the adaptation region established in block 41 . As long as these filtered parameters 
remain within the adaptation region, nothing is done. If either of the filtered parameters fall outside the adap- 
40 tation region, however, a new effective burst length is computed in block 43 and the new effective burst length 
is used to request a new connection for this user's traffic. Before proceeding to a more detailed description of 
the adaptation process, the following variables wID be defined: 

R The maximum bit rate, in bits per second, of the input traffic asrequested by the user source to initiate 
the connection. 

45 m The mean bft rate, in bits per second, of the input traffic as requested by the user source to initiate the 
connection. 

b The mean burst length, in bits, of the input traffic as requested by the user source to initiate the con- 
nection. 

t The sampling period of both of filters 45 and 47. The filters 45 and 47 receive measurements and report 
so filtered outputs to the compare unit 42 every t seconds, for a succession 1 , 2,- n, - periods, at which 

time compare unit 42 makes a decision. 
m„ The raw measurement of the mean bit rate of the input traffic for the nth sampling period of duration t 
\n The raw measurement of the red marking probability being used in the leaky bucket 34 during the nth 
sampling period of duration f. 
55 m n The filtered value of the mean bit rate, as filtered by bit rate filter 45 of FIG. 4, for the input traffic at 
the end of the nth sampling period. 

The filtered value of the red marking probability, as filtered by red marking probability filter 47 of FIG. 
4, for the leaky bucket at the end of the nth sampling period. 



5 



EP 0 643 514 A2 



b% The effective burst length of the Incoming traffic at the end of the nth sampling period, used to request 
a new connection. 

-f The green token generation rata used in the leaky bucket module 34 of FIG. 3 during the nth sampling 
period. The green token rate determines the rate at which packets marked green can be injected into 
5 the network. 

/WJ The size of of the green token pool in the leaky bucket module 34 of FIG. 3 during the nth sampling 
period. The size of the green token pool determines the length of green packets injected into the net- 
work. 

Measuring the mean bit rate m n of the incoming traffic is simple. A counter counts the number of bits re- 
10 cetved during the sampling period t and divides this number by the length t. Similarly, the red marking prob- 
ability ^„ is equal to the number of packets marked red during the sampling period t divided by the total number 
of packets transmitted during the period f . It is these raw figures which are delivered to filters 45 and 47, re- 
spectively, every t seconds. It will be noted that adaptation region (established when the connection is initially 
set up and once per bandwidth adjustment thereafter) is established in terms of the units of the raw measure- 
is ments m n and to avoid unnecessary calculations. Similarly, the effective burst length b# is also calculated 
only if the filtered A fl or £ n falls outside of the adaptation region^ As will be discussed hereinafter, an adaptation 
region is established around a target red marking probability % T and a target mean bit rate m T corresponding 
to the mean bit rate for the sampling interval / after which the previous update was implemented and which 
has been used by the network since that previous update. 
20 As taught in the aforementioned patent application Serial Number 942,873, the effective burst length can 
be calculated by substituting for the actual traffic stream an equivalent traffic stream having the same peak 
rate R, the same mean bit rate m, and the same mean burst length b, but where the substituted traffic stream 
conforms to a model having an on/off process where the on and off periods are independent (mutually and 
among themselves) and are exponentially distributed. This "exponential substitution* process is used to cal- 
25 culate an effective burst length for the actual traffic stream so that the actual traffic steam has the same packet 
loss probability e if the substituted traffic stream were fed to the same transmission link. The packet loss prob- 
ability e for the substituted traffic, as taught in the afore-mentioned application, is given by 
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_ R(c - m)exp(-n//) ttf ) 

(R - m)c - m(R - c)exp(- tUv) k } 



where 

X(c - m)FP 



(2) 



1 (R - m)(R - c)c 

where c is the speed of the transmission facility being used and X is the buffo' size of the transmission facility. 
35 Solving equations (1) and (2) for effective burst length gives 

»,-i|nn *<«-») + y-») <3) 

8C(K - m ) 

In the adaptation system of FIG. 4, the transmission facility is the leaky bucket module 34 and hence the 
packet loss probability e is the red marking probability ^ and the buffer size X is the size of the green token 
pool M n . That is, the effective burst length is given by 

5»if(R - m fl ) 

where 

„ _ AW - A,)* 2 

(R - m n y(R - ■fyf' (5) 

If the loss (red packet} probability fails outside of the desired range, it is together with R and m m which 
are passed to the connection agent 38 of FIG. 3 in order to attempt a bandwidth update for the connection. If 
the bandwidth update request is accepted, new leaky bucket parameters y"* 1 and Af * 1 are computed based 
on the new traffic parameters R, m n and b* The calculation and use of these leaky bucket parameters is taught 
in detail in the aforementioned copending application Serial Number 07/943,097. 

The filters 45 and 47 report estimates of mean bit rate and red marking probabilities, respectively, every 
t second. The value of t is determined by how fast connection agent 32 can get replies to a request launched 
55 on the network since it is useless for filters 45 and 47 to supply new estimates which might require new con- 
nection requests faster than such requests can be processed. Thus the sampling rate t is dependent on the 
maximum round trip delay of the network 10 of FIG. 1 and thus dependent on the network implementation. 
Each of filters 45 and 47 maps the current raw measurement and all of the previous measurements into an 
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estimate of the filtered value. Let x 1v x* x„ be the raw measurements and x t , x 2 , .... x„ be the estimates 
(where x is either /n or £). While the mapping of filters 45 and 46 can be any function, in the preferred embodi- 
ment of the present invention this mapping is exponential. That is, the nth estimate x„ is given by 

x„ = ox„. 1 + (1 - a]x„. (6) 

5 where the filter parameter a, lying between zero and one (tXa< 1 ), determines the relative reliability of the two 
terms in equation (6). The value of a is determined as follows, taking up the mean rate filter 45 of FIG. 4 first 
Associated with the mean bft rata m is a number T m which is the amount of time required to collect the smallest 
amount of information from which a "sufficiently accurate" estimate of m can be made. Let N m be the number 
of raw measurements of m received by filter 45 in time T m i.e., N m is the smallest integer larger than TJt N m 

10 then is the smallest number of raw measurements of m needed to determine a m Nm which is statistically reliable. 
For this reason, when it is first initialized, filter 45 does not report any estimates to compare unit 42 until N m 
measurements have been received. N t and T K can be defined similarly to represent the minimum number and 
minimum collection period of raw measurements for statistically reliable red marking probabilities in filter 47. 
If all of the mean rate estimates are constant and equal to m, the N m th estimate is given by 

15 A W( .= <^ + (1 - <#)<m - fa. (7) 

If, for example, it is desired to keep m Nm within 90% of m from the initial condition mso, then a„ satisfies 

1 - c£ = 0.9 <8) 

20 and 

a* = ai«V (9) 

Using confidence interval analysis, the determination of the value of T m (and of 7"J proceeds in the following 
steps: 

1. Identify independent (or almost independent) identically distributed (II Ds) "experiments" involved in 
25 measuring mean rate; 

2. Determine the minimum number of experiments needed to achieve a desired confidence interval; and 

3. Determine T m as the amount of time needed to collect the minimum number (from the previous step) 
of experiment results. 

Assuming that the traffic can be modeled as an on/off process satisfying the independence and exponerv 
30 tial distribution assumptions, then the on/off cycles are the IID experiments. Let B n and i n be the nth on and 
off times, respectively, and let the mean bit rate in the nth cycle of length B D + /„ be Y n = RBfj(B n + /„). The 
means of the common exponential distributions of B n and /„ are represented by p*~ 1 and nT 1 , respectively, and 
the mean of the IID random sequence {Y m n ■ 1, 2, •»} is given by 
^ rrij = pR where p ■ n//(n, + \L b ) (10a) 

and its standard deviation is approximately given by 

a, = V2p(1 - p)R. (10b) 

Let M m be the smallest sample size in the sequence Y n such that the sample mean /nhasa confidence interval 
(m-Zm, m+Zm) and a confidence level 8 m . The let w m be obtained from the normal distribution table to satisfy 
40 P{-w m £W*w m } = 0 mt (11) 

where IV is normally distributed with zero mean and unit variance. Then M m is given by 

W m = (^^)2. (12) 

45 The time interval T m required for obtaining the desired confidence level is then as the time required to obtain 
M m raw measurements. If Is the mean burst length, then one measurement takes, on the average, bjm 
seconds and T m is given by 

T m = M m ^. (13) 

60 For the red marking probability filter 47 of FIG. 4, it i3 assumed that every packet has the same probability 
of being marked red and that whether or not it is marked red is independent of all of the others. Although this 
is not strictly true, particularly when multiple packets arrive back-to-back during an on period, the results used 
herein are not very sensitive to the confidence level of the red marking probability £ estimates. Let X, be a 
random variable indicating whether the Hh packet is marked red (X/ = 1) or green (X) = 0). Then 
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y — I 1 with probability £ 

\ 0 with probability (1 - 0 ' < 14 ) 

5 

The standard deviation a, of X, is 

<rx = V$ - p. (15) 

Using the same steps used to compute the mean bit rate measuring period T m , the smallest number A^ofX 
10 observations to achieve a confidence interval of - z^ t % + z§) and confidence level 0$ is given by 

A* = <^) 2 . (16) 

where is obtained similarly to The desired measuring period T K is then obtained as the time to observe 
the M$ packets. If L is the mean packet length (in bits), a packet will arrive on the average every L/m seconds 
15 and T$ is given by 

T t = M t £. (17) 

The filters 45 and 47 of FIG. 4 thus utilize the measurement periods T m and T K defined in equations (1 3) 
and (17), respectively. Filtered values A n and ^ of the mean bit rate and the red marking probability thus de- 
livered to compare unit 42 once every T m and T K seconds, respectively. These filtered values form a set which 
can be compared in unit 42 to the acceptable sets of these values to determine whether or not adaptation is 
required, i.e., a new connection request is warranted. This comparison will be discussed in connection with 
the adaptation region disclosed graphically in FIG. 5. 

Before proceeding to a discussion of FIG. 5, it should generally be noted that a new connection is requested 
when the mean bit rate or the "burstiness" of the input traffic renders the leaky bucket mechanism 34 incapable 
of insuring 1) the negotiated overall packet loss probability of that connection or 2) the packet loss probability 
of other connections sharing the same links. On the other hand, a new connection can also be requested when 
the mean bit rate or the burstiness of the in put traffic falls sufficiently low that 3) the throughput of the network 
can be increased by reducing the bandwidth allocation for that connection, in the process also reducing the 
user cost based on allocated bandwidth. In accordance with the present invention, these three conditions are 
used to define the bounds of the adaptation region shown in FIG. 5, as will be shown below. 

In FIG. 5 there is shown a graphical representation of the adaptation region in accordance with the present 
invention for adapting the bandwidth of a packet network connection in response to changes in the mean bit 
rate and/or the burstiness of the traffic input to the network. Curve 51 , between points 50 and 52, represents 
the limit on the current connection's overall packet loss probability. For the sake of simplicity, if the overall 
packet loss probability % of this connection stays within the negotiated order of magnitude, that is, if e c e (0.5, 
5.0) x e n then no adaptation adjustment is necessary, if Sg and «, are the loss probabilities, respectively, for 
the green and the red packets in the network, it is assumed that the bandwidth allocation is done to equalize 
these two loss probabilities and that the buffer management satisfies the order relationship 

0(«J* = 5rO(e r ) + (1 - trJOfe^ (18) 

where O(e) is defined as 1 0- n when (0.5,5)*1 Q- n where (n is an integer) 

and 4 r is the red marking probability, i.e., $ T = O («g yo (e, ). When % increases beyond £ r , the overall 
loss probability becomes unacceptable and the bandwidth must be adjusted upward. In FIG. 5, curve 51 is 
the collection of (/n, £>J points that yield % = k& T where is 5, the high order of magnitude limit on the target 
probability loss. That is, curve 51 is the collection of sets of values of m and to* which satisfy 
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SO 



(R - m)y - m(R - y)exp(- t^d*) 

where 

M a (y - m)& 



-*i5* (19) 
(20) 



1 (R - mX* - 7)f" 

Curve 53 between points 52 and 54 of FIG. 5 represents the I imit on the amount this connection is allowed 
to unfairly impact on other connections sharing the same transmission links. It is first noted that launching of 
55 excess traffic into the network is more likely to first affect the average delay seen by other connections before 
impacting the loss probability of these other connections, due mainly to the red discard thresholds imposed 
at the intermediate nodes. The buffer occupancy of the buffer in a link is a measure of this average connection 
delay when this connection source is Input into the buffer with an effective transmission rate equal to the band- 
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width reserved on that link for this connection. The presence of the leaky bucket is ignored since, if curve 51 
is satisfied, then all of the traffic will enter the network without being shaped by the leaky bucket, as long as 
there is no limit on the amount of red traffic that can be sent into the network. The buffer content due to both 
the green and the red traffic must be considered since the traffic priority is on a first-come, first-served, basis 

5 and the impact of the red discard threshold on the average delay will be minimal since the red discard threshold 
is designed to control discard probabilities. 

Curve 53 of FIG. 5 is determined in the following way. Let v be the average buffer content for the original 
source. When the average buffer occupancy increases to a multiple k 2 xv(k 2 > 1), the traffic is deemed to be 
unacceptable to the network. The value of the coefficient k 2 depends on the bandwidth of the connection rel- 

10 ative to the link's capacity, the larger this ratio, the smaller the value of Moreover, since a connection path 
generally consists of more than one link, the buffer capacity in some links may become unacceptable sooner 
than the others. The link in which the buffer occupancy first becomes unacceptable is that link in which the 
connection's reserve capacity contribution is the smallest As discussed in the aforementioned copending ap- 
plication Serial Number 943,097, the green token generation rate y e of the leaky bucket is set equal to this mini- 

15 mum reserve capacity contribution. Therefore, curve 53 is formed by the collection of all (/n, bj) points which 
satisfy the relationship 



20 



25 



40 



where 



g = * 2 v, (21) 



5 = <Y ~ . (22) 



Curve 56 between points 55 and 57 represents the lower limit on the adaptation region of FIG. 5 and is 
heavily dependent on the customer billing structure. For each connection, a different amount of bandwidth Is 
reserved on each link along the path, depending on the initially assumed characteristics of the traffic. The 
amou nt charged to the customer is, in general, related to these amounts of reserved bandwidth. If it is assumed 
that the customer is billed by an amount proportional to the bandwidth y, then curve 56 is the collection of all 
(m, bj) points that result in the y parameter being some fraction /fe of the current value, where the fraction k$ 
depends on the size of the connection and implementation constraints such as bandwidth update overhead. 
30 For lower bandwidth connections, /r 3 can be smaller since update overhead is the same regardless of the band- 
width of the connection. (???) 

For simplicity, it is assumed the all of the buffer sizes are the same (X) and that the target loss probabilities 
for each link are the same (*). It is further assumed that the traffic characteristics are such that y is equal to 
the equivalent capacity of the connection. In this case, curve 56 of FIG. 5 is the collection of ail (m, bj) points 
35 satisfying the relationship 

where y is the current token generation rats and 



r . ffl^Ei^lnl. (24) 



In accordance with the illustrative embodiment of the present invention, the adaptation region 56 of FIG. 5 is 
bounded by equations (1 9)-(24) and represents a region within which adaptive adjustments in the connection 
bandwidth are not required. Above region 58, the bandwidth is adjusted upward while below region 56 the band- 
45 width is adjusted downward. The values of k u k* k s can be adjusted to reflect local requirements of the network 
1 0 of FIG. 1 . Generally, however, fci lies in a range between 2 and 1 0, k 2 lies in a range between 1 .1 and infinity, 
and k$ lies in a range between zero and 09. In one application, ^ « 5, k 2 ° 2 and k$ = 0.8 performed well. 

It will be noted that the comparisons necessary to determine violations of the adaptation requirements 
summarized in FIG. 5 require conversion of the red marking probability data into effective burst lengths to per- 
se mit the comparison. Since this conversion is a complex floating point operation and is required every compar- 
ison period, considerable time would be lost in calculating effective burst lengths which would neverbe used 
(since this effective burst length is otherwise used only if a new connection is requested, i.e.. the m m ^ points 
foil outside of the region 58 of FIG. 5). In accordance with the present invention, the computation load on the 
bandwidth management subsystem of the present Invention is significantly reduced by remapping the adap- 
55 tation region 58 of FIG. 5 in terms of the actual measured quantities m„ and This remapping is shown in 
FIG. 6. 

In FIG. 6 there is shown a graphical representation of the adaptation region of FIG. 5 mapped to the average 
bit rate (m) and red packet marking probability {Q plane. The same reference numerals used in FIG. 5, with 
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an accent attached, are used in RG. 6 to represent corresponding portions of the adaptation region. Thus adap- 
tation region 58' of FIG. 6 is bounded by curves 51', 53' and 58', corresponding to equations (19), (21) and 
(23), respectively. Any point on the (/n, bj) plane can be mapped into a point on the (m. Q plane, given the 
initial traffic characterizations R, m, bj) using the relationship 
s R<y - m)exp(-n/^) 

where r\ is defined in equation (20). Although this remapping also involves complex computations, it need be 
undertaken only once every bandwidth update rather than once every decision period (r). Hie set of points 
within region 58' can be denoted by R(m ; , % T ), where A, is the m value used by the subsystem of FIG. 3 to set 
10 up the currently existing connection and £ r is the target red marking probability currently being enforced by 
the leaky bucket module 34. 

In RG. 7 there is shown a flow chart of the operation of the estimation and adaptation module 33 of FIG. 
3, shown in more detail in RG. 4. The adaptation algorithm shown in FIG. 7 can be summarized as follows: 

15 ACCESS CONTROL ALGORITHM 



A Initialization: Given the initial traffic parameters (R t i>£); 
1. Compute N m and ou using equations (9)-(13). 
20 2. Compute N$ and c^, using equations (15)-(17). 

3. Compute the region R(m 0 , in the (m, £) plane. 

4. Set m 0 = mo, $o = 5r. " = 1 andy = 0. 

B. At the end of each measurement interval n, n - 1 , 2, •»-, filter the measurements of m and using the 
exponential filter of equation (6). 
25 C. Test for (n < A/ m ), i.e. has a sufficient number of raw m measurements been received to calculate a suf- 
ficiently accurate m„? 

1. If yes, test for (n<j+ i.e., has a sufficient number of raw £ measurements been received since 
the last connection establishment to calculate a sufficiently accurate 

a. If not, increment n = n + 1 and return to step (B). 
30 b. If yes, test for (rhp ^) € R (m /f !■/), i.e., is the current (m t Q set within region 58' of FIG. 8? 

• If within region 58', increment n and return to step (B). 

• If outside of region 58\ skip to step (D), requesting a new connection, using the parameters 

2. If (n £ Nm), test for (n <y + N0 9 Le. , if a sufficient number of raw m measurements has been received, 
35 has a sufficient number of raw £ measurements been received since the last connection establishment? 

a. If not, increment n and return to step (B). 

b. If yes, test for (m„, £ n ) e R (m> £7). i-e., is the current m, % pair within region 58' of RG. 6? 

• If within region 58', increment n and return to step (B). A 

• If outside region 58*, skip to step (D), requesting a new connection, using the parameters (m* y . 
^ D. Request connection agent 38 (FIG. 3) for a new connection, using the parameters received from one 

of the above steps. 

E. Wait for new system parameters (y, M g , Not, N fy c^) from connection agent 22. This wait is required 
only for bandwidth increases. Bandwidth decreases can always be accommodated within the originally re- 
served bandwidth. 
^ F. Reset all system parameters to the new values. 

1. Set = £7 and set/ = n. 

2. Increment n a n + 1. 

3. Return to step (2). 

Referring to FIG. 7, starting at start box 70, box 71 is entered where the connection parameters from the 
50 Initial connection request for a connection (R, m 0 , and b* are received (at the same time the initial connection 
request message is launched on the network by connection agent 32 (FIG. 3)). In box 72, all of the variable 
parameters of the estimation and adaptation module of FIG. 4 are computed. In particular, the minimum number 
N m of raw mean bit rate measurements, the mean bit rate filter 45 weighting factor o^, the minimum number 
A/ ( of raw red token probability measurements, and the red token probability filter 47 weighting factor a t . In 
55 addition, the initial filtered value of m 0 is set at the initial request value of m 0l the initial filtered value of £0 is sat 
at the target red packet marking probability £ r in box 72. Finally, in box 72, the measurement period index n 
is set to "1" and the connection request number j is sat to "0". In box 73, the next raw measurements are re- 
ceived, a mean bit rate m A and a red token probability These raw values are used In box 74 to calculate the 
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filtered values m„ and using the exponential filter equation (6). 

The output of box 74 is applied to decision box 75 to determine if a sufficient number of mean bit rate raw 
measurements have been received to guarantee a sufficiently accurate mean bit rate estimate. If not, decision 
box 76 is entered to determine if a sufficient number of red marking probability raw measurements have been 

s received to guarantee a sufficiently accurate red token marking probability estimate since the last connection 
update. Note that the raw mean bit rate measurements need be monitored for a minimum number only once 
for a connection, wh ile the red token marking probability measurements must be mo nitored for a minimum num- 
ber after each connection update since the leaky bucket parameters are updated whenever a connection is 
updated. If neither sufficient mean bit rate measurements nor sufficient red token marking probability meas- 

w urements have been received, as determined by decision boxes 75 and 76, box 84 is entered where the meas- 
urement index n is incremented by one. Box 73 is then re-entered to receive the next raw measurement values. 

If a sufficient number n of raw red token marking probability measurements have been received (at least 
j + N£, as determined by decision box 76, decision box 78 is entered to determine if the (A^J point in the 
(/n,£) plane is within the adaptation region 58' of FIG. 6. If so, no adaptation is required and box 84 is entered 

is to increment n and re-enter box 73 to get the next raw measurement If the point (mfa) is outside of the adap- 
tation region 58' of FIG. 6, as determined by decision box 78, box 79 is entered to request a connection with 
the new traffic characteristics (/tyJU). That is, if sufficient 2; samples have been received to accurately calculate 
the red token probability estimate, but insufficient samples have been collected to accurately estimate the 
mean bit rate, nevertheless, the new ^ estimate is combined with the mean bit rate in force from the last oon- 

20 nection update and compared to the adaptation region 58* of FIG. 6. If these values fall outside of the adaptation 
region 58', then a new connection is requested in box 79. If a new connection is negotiated, a new adaptation 
region is calculated, new leaky bucket parameters are calculated, and new filter values are calculated. The 
details of this process in boxes 79, 82 and 83 will be discussed in connection with FIG. 8. After the connection 
is updated, box 84 is entered to increment n and box 73 re-entered to await the next measurement sample. 

25 If it is determined in decision box 75 that a sufficient number of raw m„ samples have been received to 
provide a sufficiently accurate estimate of the mean bit rate, decision box 77 is entered where it is again tested 
as to whether a sufficient number of raw red marking probability samples have been received since the last 
leaky bucket update to estimate the red marking probability. Even if insufficient ^ samples have been re- 
ceived, decisiop box 81 Is entered to determine if the (ro^) point falls within or outside of the adaptation region 

30 58'. If the point does fall outside of the adaptation region 58', as determined by decision box 81 , box 
83 is entered to update the connection, using the new mean bit rate m„ but the target red marking probability 
% T as the update parameters. That is, even when sufficient raw ^ samples have no? been received, neverthe- 
less the connection is updated if the mean bit rate has changed so much that the (A^y is no longer in the 
adaptation region 58'. After the connection is updated in box 83, box A 84 is entered to increment n and box 73 

35 re-entered to await the next measurement value. If the test point {m^) is within the adaptation region 58', as 
determined by decision box 81 , box 84 is entered to increment n and box 73 re-entered to await the next meas- 
urement value. 

If it is determined in decision box 77 that sufficient ^ values have been received to estimate the red mark- 
ing probability, decision box 80 is entered to test if the (m„£ T ) point is within the adaptation region 58'. rf so, 
40 box 84 is entered to increment n and re-enter box 73. If the (A„,| T ) point is outside of adaptation region 58', 
box 82 is entered to update the connection with the parameters (m>^). Box 84 is then entered to increment n 
and return to box 73. 

The procedure shown in FIG. 7 is, in the preferred embodiment. Implemented by programming a general 
purpose computer. If some of the highly iterative portions of the algorithm of FIG. 7 which are required on a 

46 per connection basis are implemented with high-speed special purpose circuitry, it is important that the inter- 
face between the high-speed circuitry and the computer be designed to min imlze computer Interrupts. It should 
also be noted that the adaptation region 58' of FIG. 6 can be implemented by sampling the value of £'s at a 
number of m's covering a range of m values to create a two dimensional array fc L {mft t % H ( m d* 1 B 1 . 2, .... K} 
of regular m values and the corresponding £ values for each of curves 51', 53', and 56*. The m values need 

bo be sampled only In the range between zero and the currently used green token source rate y since (m„ > y) 
always implies the need fa- an upward bandwidth adjustment (y is the maximum rate at which data can flow 
through the leaky bucket). The values of other £s 5a n be interpolated between the stored values. Given the 
resulting^ and £, values, the relative position of ^ with respect to and ^(rh R ) can be determined. If 
the new t, estimate lies between these two values, the new measurement values are within region 58'; other- 

55 wise the values are outside of the adaptation region 58'. 

One modification of the algorithm of FIG. 7 is to decrease the bandwidth more gradually than the bandwidth 
is increased. This can be accomplished by introducing an hysteresis parameter h (0 < h< 1) into the compu- 
tations for m„ and Aw for estimates below the lower boundary curve 56' of adaptation region 58*. That is, if the 
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estimates (m n Xd are telow curve 56\ then 

m 0 4-hmj+(\ 'h)fh„ (26) 

and 

(1 -h)$a (27) 

5 

are used in computing the request vector and the new leaky bucket parameters. The value of h can be set at 
or near 0.5. Another variation of the algorithm of FIG. 7 is to delay requests for updates for a larger number of 
raw measurement samples, providing a more conservative adaptation strategy. 

In RG. 8 there is shown a flow chart of the process for updating a connection, shown generally In boxes 
10 79, 82, and 83 of RG. 7 and in step (D) of the access control algorithm outlined above. The update connection 
algorithm of boxes 79, 82 and 83 of FIG. 7 can be summarized in the following substeps, utilizing the values 
of m and 5 received from the flow chart of FIG. 7. 

CONNECTION UPDATE ALGORITHM 

15 

A Compute A* from equation (3) and send a bandwidth update request message (FIG. 2) including the 
new bgf, m and R. 

B. Walt for a reply from the remote connection agent If this is a bandwidth reduction request, no waiting 
is required; proceed directly to step (C). 

20 1 . If the reply is negative, terminate or reroute the connection. 

2. If the reply is positive, proceed to step (C). 

C. Compute the following quantities for the new connection: 

1. Compute the leaky bucket parameters y (the green token source rate) and M g (the green token pool 
size). 

25 2. Compute the new adaptation region R( m£ T ) (RG. 6). 

3. Compute new values for W„, a*,, and a$. 

Referring more particularly to FIG. 8, starting at start box 90, box 91 is entered where the value of the new 
effective burst length /> rf is calculated, using equation (3). Box 92 is then entered to launch a request from con- 
nection agent 32 (FIG. 3), using the connection characterization values specified in box 79, 82 or 83 of FIG. 

30 7. Decision box 93 is then entered to determine if the request is for a bandwidth increase or a bandwidth de- 
crease. If a bandwidth increase is requested, box 94 is entered where the system waits until the connection 
request propagates along the connection path to the remote receiving agent and a confirmation of the con- 
nection is returned to connection agent 32. Decision box 95 is then entered to determine if the reply is positive 
(update granted) or negative (update denied), if the update is denied, box 96 is entered to terminate the con- 

35 nection or to invoke other facilities, not herein disclosed, to reroute the connection. The procedure terminates 
in box 97. 

if the connection request reply is positive, as determined by decision box 95, box 98 is entered where new 
leaky bucket parameters y (green token source rate) and M g (green token buffer size), as taught in the copend- 
ing application Serial Number 943,097, identified above. If a bandwidth decrease is requested, as determined 

40 by decision box 93, box 98 is entered directly since no wait for a response to the connection request message 
is required; the bandwidth will always be available from the wider bandwidth connection being replaced. After 
calculating the new leaky bucket parameters in box 98, box 99 is entered to compute the new adaptation region 
58'. using equations (19). (21) and (23), and, using equation (25) translated to the m& plane. In box 100, new 
filter values a mt /^.a^are computed which control filters 45 and 47 during the upcoming connection update. 

45 It can be seen that the access control mechanism of the present invention insures that adequate bandwidth 
will be provided for each connection even when the characteristics of the signal source associated with that 
connection change significantly. This "dynamic adaptation" to the characteristics of the signal source, however, 
is controlled to prevent overly rapid updates to the connection parameters and the network overhead associ- 
ated therewith. The bandwidth, measurement filters and leaky bucket access mechanism are all changed to- 

60 gether to insure compatible treatment of the Input traffic stream. 

One obvious modification of the present invention is to replace the red marking probability £ of the leaky 
bucket mechanism with the probability that the token level at the leaky bucket mechanism is less than some 
predetermined threshold level. This latter metric is useful when a very small red marking probability is desired, 
for example if the red loss probability In the network were very high. 

55 
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Claims 

1. A packet communications network for interconnecting source nodes and destination nodes for establish- 
ing a connection for the transmission of digital packets of information from a source node to a destination 
5 node in said network, said communications system including a dynamic access control mechanism com- 

prising 

means for monitoring the mean bit rate of signals from said source, 

a leaky bucket control circuit for controlling the flow of said signals from said source into said net- 
work, 

10 means for monitoring the loss probability of packets introduced by said leaky bucket control circuit, 

means for establishing boundaries of corresponding pairs of measurements from said mean bit rate 

monitoring means and said loss probability measuring means, and 

means responsive to a pair of said mean bit rate and loss probabii ity measurements falling outside 

of said boundaries for updating the bandwidth of said connection. 
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2. The packet communications network according to claim 1 wherein said means for establishing boundaries 
comprises 

means for determining an order of magnitude range around said mean bit rate and said loss prob- 
abilfty measurements. 

3. The packet communications network according to claim 1 wherein said means for monitoring said mean 
bit rate comprises 

filtering means for filtering a plurality of said mean bit rate measurements. 

4. The packet communications network according to claim 1 wherein said means for monitoring said loss 
probability comprises 

filtering means for filtering a plurality of said loss probability measurements. 

5. The packet communications network according to claim 2 wherein said boundaries represent values of 
the mean bit rate m and the packet loss probability % satif ying the relationships 

R(Y - m)exp(-n/*W) - k E 

(R - m)y - m(R - y^pC-n/i)*) 1W 

where 
and 

f 

where 

5 = (y-m)R* 

and 



45 where y is the current token generation rate and 

R is the maximum bit rate of signals from said source node, 
7 is the green token source rate of said leaky bucket control circuit, 
so 6^ is the effective burst length of signals from said source node, 

k\ is a constant between two and ten, 

%r is the target red token loss probability of said leaky bucket control circuit 
M g is the size of the green token buffer in said leaky bucket control circuit, 
k 2 is a constant between 1.1 and infinity, 
55 v is the average buffer content of signals at said source node, 

a 7 is the target loss probability, 
* 3 is a constant between zero and 0.9. 
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6. A method of dynamically adapting access to a packet communications network for interconnecting source 
nodes and destination nodes for the transmission of digital packets of information from said source node 
to said destination node, said method including the steps of 

monitoring the mean hit rate of signals from said source node, 
5 uti llzing a leaky bucket control circuit to control the flow of said signals from said source node into 

said network, 

monitoring the loss probability of packets introduced into said network by said leaky bucket control 

circuit, 

establishing boundaries of corresponding pairs of measurements from said mean bit rate monitor- 
10 ing step and said loss probability monitoring step, and 

means responsive to a pair of said mean bit rate and loss probability measurements felling outside 
of said boundaries for updating the bandwidth of a connection between said source node and said des- 
tination node. 

16 7. The method according to claim 6 wherein said step of establishing boundaries comprises the step of 

determining an order of magnitude range around said mean bit rate and said loss probability meas- 
urements. 
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8. The method according to claim 6 wherein said step of monitoring said mean bit rate comprises comprises 
the step of 

filtering a plurality of said mean bit rate measurements. 

9. The method according to claim 6 wnerein said step of monitoring said loss probability comprises com- 
prises the step of 

filtering a plurality of said loss probability measurements. 

10. The method according to claim 7 wherein said step of establishing boundaries further include the step of 

determining values of said mean bit rate m and said packet loss probability % satifying the relation- 
ships 

R(y- /n)exp(-Ti/b„) u 

(R - m)y - m(R - T )axp( w 



where 



and 



where 



and 



(R- m)(R- y yf' 



(y- myp 

(R- mXR-yhbJ 



(22) 



where y is the currant token generation rate and 

r-^^WlrA (24) 

R is the maximum bit rate of signals from said source node, 

7 is the green token source rate of said leaky bucket control circuit 

b rf is the effective burst length of signals from said source node, 

so ki 'isH constant between two and ten, 

5r is the target red token loss probability of said leaky bucket control circuit, 
M g is the size of the green token buffer in said leaky bucket control circuit, 
k 2 is a constant between 1.1 and infinity, 
v is the average buffer content of signals at said source node, 

55 8 T is the target loss probability, 

k$ is a constant between zero and 0.9. 
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FIG. 3 
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